clc  
clear all
close all
%% input
 %mydata = load(y10)
% load y10
% Case=y10;
tic
cd ('E:')
load Betaband_d212posteo
m=3;
tau=1;
X=Betaband';
mydata=Betaband';
%X=Case;
% 3D State Space of EEG signal   
% [m,tau]=optimum_m_tau(X);%compute optimum dimension and tau
[X1,T1]=phasespace(mydata(:,1),m,tau);  %X1 indicate channel1
%[m,tau]=optimum_m_tau(X)
[X2,T2]=phasespace(mydata(:,2),m,tau);
[X3,T3]=phasespace(mydata(:,3),m,tau);
[X4,T4]=phasespace(mydata(:,4),m,tau);
[X5,T5]=phasespace(mydata(:,5),m,tau);
[X6,T6]=phasespace(mydata(:,6),m,tau);
[X7,T7]=phasespace(mydata(:,7),m,tau);
[X8,T8]=phasespace(mydata(:,8),m,tau);
[X9,T9]=phasespace(mydata(:,9),m,tau);
[X10,T10]=phasespace(mydata(:,10),m,tau);
[X11,T11]=phasespace(mydata(:,11),m,tau);
[X12,T12]=phasespace(mydata(:,12),m,tau);
[X13,T13]=phasespace(mydata(:,13),m,tau);
[X14,T14]=phasespace(mydata(:,14),m,tau);
[X15,T15]=phasespace(mydata(:,15),m,tau);
[X16,T16]=phasespace(mydata(:,16),m,tau);
[X17,T17]=phasespace(mydata(:,17),m,tau);
[X18,T18]=phasespace(mydata(:,18),m,tau);
[X19,T19]=phasespace(mydata(:,19),m,tau);
save mydata % phasespace save

%% simulation

%SL=codeSL;
FSL=codeFSL;


% ChannelsSL1=find((SL(1)-0.005)<=SL & SL<=SL(1)+0.005);
% ChannelsSL2=find((SL(2)-0.005)<=SL & SL<=SL(2)+0.005) ;
% ChannelsSL3=find((SL(3)-0.005)<=SL & SL<=SL(3)+0.005) ;
% ChannelsSL4=find((SL(4)-0.005)<=SL & SL<=SL(4)+0.005) ;
% ChannelsSL5=find((SL(5)-0.005)<=SL & SL<=SL(5)+0.005) ;
% ChannelsSL6=find((SL(6)-0.005)<=SL & SL<=SL(6)+0.005) ;
% ChannelsSL7=find((SL(7)-0.005)<=SL & SL<=SL(7)+0.005) ;
% ChannelsSL8=find((SL(8)-0.005)<=SL & SL<=SL(8)+0.005) ;
% ChannelsSL9=find((SL(9)-0.005)<=SL & SL<=SL(9)+0.005) ;
% ChannelsSL10=find((SL(10)-0.005)<=SL & SL<=SL(10)+0.005) ;
% ChannelsSL11=find((SL(11)-0.005)<=SL & SL<=SL(11)+0.005) ;
% ChannelsSL12=find((SL(12)-0.005)<=SL & SL<=SL(12)+0.005) ;
% ChannelsSL13=find((SL(13)-0.005)<=SL & SL<=SL(13)+0.005) ;
FSL=FSL';
save FSL_Betaband_d212posteo FSL
Mean_FSL_Betaband_d212posteo=mean(FSL);
save Mean_FSL_Betaband_d212posteo Mean_FSL_Betaband_d212posteo
toc